Superkicks in Hyperbolic Encounters of Binary Black Holes 
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Generic inspirals and mergers of binary black holes produce beamed emission of gravitational 
radiation that can lead to a gravitational recoil or kick of the final black hole. The kick velocity 
depends on the mass ratio and spins of the binary as well as on the dynamics of the binary configura- 
tion. Studies have focused so far on the most astrophysically relevant configuration of quasi-circular 
inspirals, for which kicks as large as ~ 3,300 km s _1 have been found. We present the first study 
of gravitational recoil in hyperbolic encounters. Contrary to quasi-circular configurations, in which 
the beamed radiation tends to average during the inspiral, radiation from hyperbolic encounters is 
plunge dominated, resulting in an enhancement of preferential beaming. As a consequence, it is 
possible to achieve kick velocities as large as 10, OOOkms -1 . 
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Numerical relativity estimates of the gravitational re- 
coil or kick inflicted on the final black hole (BH) from 
generic inspirals and mergers of binary black holes (BBH) 
have triggered tremendous excitement in astrophysics. 
This is mainly due to the fact that most galaxies host 
a supermassive black hole (SMBH) at their centers [l], Q • 
As galaxies merge, a kick to the final BH from the coa- 
lescence of the BHs at the galactic cores could have pro- 
found implications in subsequent mergers, affecting the 
growth of SMBHs via mergers as well as the population of 
galaxies containing SMBHs. In addition, there have been 
several suggestions of direct observational signatures of 
putative BH recoils 0, i, H, H 0, H, with one stud y Hi 
presenting evidence for the first candidate of a recoiling 
SMBH. 

BH kick velocities depend on the mass ratio and spins 
of the merging BHs as well as on the initial configura- 
tion and subsequent dynamics of the binary. Studies 
published to date have concerned the most astrophysi- 
cally relevant configuration, that of quasi-circular inspi- 
rals [10, [U dM 03 ■ O ne remarkable discovery has been 
kicks of a few thousand km s~ found in configurations 
of equal-mass binaries with initially anti-aligned spins in 
the orbital plane [14|, LLa LL6|. For near extremal spins 
{a/m — 0.925), recoils as large as 3, 300 km s~ have 
been computed |17| . 

Motivated by our previous study [18|] of the final spin 
of BHs from scattering mergers of BBHs, we present 
the first extension of gravitational recoil to hyperbolic 
encounters. There are crucial differences between hy- 
perbolic and quasi-circular configurations that affect the 
kick to the final BH. For quasi-circular orbits, the emit- 
ted radiation is asymmetrically beamed, and carries lin- 
ear momentum with it. But it tends to average out 
during the inspiral [l9(, producing a modest wobbling 
and drift of the center of mass of the binary. The final 
kick arises because as the binary approaches the plunge, 



the averaging loses its effectiveness, leading to a gradual 
recoil build-up. Both numerical simulations and post- 
Newtonian studies [Hj, [2(| have confirmed the gradual 
kick accumulation during the inspiral, and, in addition, 
the studies have shown that most of the recoil is gener- 
ated during the plunge. In some instances, during the 
plunge and ring-down there is also aperiod of anti-kick 
before reaching the final kick value [2l|, [23] ■ 

The main motivation for the present study was to 
consider plunge-dominated configurations to investigate 
whether kicks comparable to those for quasi-circular in- 
spirals can be found. We were surprised to find that 
kicks larger than 3000 km s _1 are in general produced for 
spin configurations equivalent to those studied in quasi- 
circular inspirals. We show below a configuration giving 
kick velocities as large as 10,000kms _1 . Two qualita- 
tive features of hyperbolic encounters contribute to these 
larger kicks. Not only are hyperbolic encounters plunge 
dominated, but the nature of the plunge is such that it 
enhances the beaming of radiated linear momentum [23j . 

We use the same computational infrastructure and 
methodology as in previous studies [1JJ, [l2|, H3 > namely 
a BSSN code with moving puncture gauge conditions. 
The hyperbolic encounters are initiated with Bowcn-York 
initial data [2J]. The data consist of two equal-mass 
BHs with masses m = AI/2 located along the x-axis: 
BH± is located at x = ±5 M and has linear momentum 
P± = ±(P cos 0, P sin 9, 0) with 8 the angle in the orbital 
plane with respect to the a;-axis. The total initial orbital 
angular momentum is then L/M 2 = 10 (P/M) sm6 z. 
The spins of each BH are in the orbital plane: BH± 
has spin S± = ±(5 cosc/)±, S sin^>±,0), with <f>± the 
angle in the orbital plane with respect to the ±x-axis. 
The parameter space of our encounters is quite large: 
{P, 9, S, <f>±}. However, from exploratory runs we have 
gained a good understanding of the parameter space and 
isolated those parameters that can be kept fixed without 
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TABLE I: Configuration parameters P/M, <j)+ and initial or- 
bital angular momentum L/M 2 for the cases with cj>- = 0° 
and S/M 2 = 0.2. The last four columns show the correspond- 
ing kick velocities of the final BH in km s _1 as well as the % 
of energy radiated. 



seriously compromising the goals of the study. 

For most of the runs, we have kept the spin magnitudes 
at S/M 2 = 0.2 (a/m = 0.8), with the exception of those 
runs used to investigate the dependence of the kick on S. 
We kept fixed also the impact angle at 6 = 153.4°. We 
considered some other angles but found that this is the 
angle for which we obtained the largest kicks. Finally, 
for most cases we kept the spin direction of BH_ located 
at x = —5 M fixed at 0_ = 0° or 45°. Once other pa- 
rameters were fixed from among a small set of values, the 
parameter we varied in general was the linear momentum 
magnitude P. 

Kicks are computed from a surface integral [2a, [261] in- 
volving the Weyl curvature tensor ^4. Although strictly 
speaking this kick formula must be evaluated in the 
limit r — > 00, we applied it at extraction radii r/M = 
{40,50,75,85,100}. The resultant kicks were fitted to 
) , &-„ i„ ^a \rt„\ - yV) + Kl j r + 



both V(r) = V£ J + K /r and V(r) 
K^/r 2 . The extrapolated r — > 00 kicks and their er- 
rors were estimated from Voo = (Vob + Vob )/2 and 
SVoc = Voo — Vob I, respectively. Unless explicitly noted, 
all the reported kick velocities and energy radiated were 
obtained from simulations with resolution h = M/0.8 on 
the mesh used for the kick computation. Every run had 
10 levels of factor-of-2 refinement, with outer boundaries 
at ~ 320 M. We discuss below the convergence of kick 
estimates as the grid spacing is decreased and the extrac- 
tion radius is increased. The estimated errors in all the 
kick results presented here due to the finite differencing 
grid spacing h = M/0.8 and extrapolation r — > 00 are of 
the order of a few hundred kms -1 . 
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FIG. 1: Magnitude of the kick velocity as a function of the 
magnitude of the initial BH spin (a/m) for the case (P/M = 
0.2665, <h+ =45°) 



Table U gives the components of the recoil along the 
initial orbital angular momentum V\\ and in the initial 
orbital plane V± in km s~ (and the total recoil V) for 
the cases with 0_ =0°. Notice that the largest recoil in 
this case happens when 0+ = 45°. Another important 
observation is that, although the dominant component 
of the kick is V\\ , as we increase 0+ a substantial compo- 
nent of the kick is also generated in the orbital plane (see 
V± ) . The rightmost column in Tabic |T] gives the r — > 00 
extrapolated radiated energy as a percentage of the to- 
tal initial energy of the binary. We see that large angular 
momentum increases the energy radiated, up to very sub- 
stantial values, as large as 13%. For even larger angular 
momenta, we expect a falloff of the radiated energy. 

For the case P/M = 0.2665,0+ = 45° in Tabic [fl we 
have carried out simulations to investigate the depen- 
dence of the kick with the initial BH spins (a/m) . The 
results are displayed in Fig. Q] We find that, to first or- 
der, the kick is proportional to a/M^, as found in quasi- 
circular orbits [13]. However, as the initial spin of the 
BHs grows, we found hints of the quadratic spin depen- 
dence also obtained in quasi-circular orbits [271 ]. 

The configuration we have found to yield the largest 
kick has the spins anti-aligned <j>+ = <j> = <p, as in the 
case of quasi-circular orbits [14|, [l5|, [l6| , and initial lin- 
ear momentum P/M = 0.3075, corresponding to initial 
angular momentum L/M 2 = 1.275. The kicks are es- 
sentially along the direction of the initial orbital angular 
momentum. Figurc[5]shows the kick velocity Vji as a func- 
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FIG. 2: Kick velocity Vji as a function of "& — — 4> for the 
cases with spin orientation <f>+ = cj>- = <f> and initial linear 
momentum P/M = 0.3075, corresponding to initial angular 
momentum L/M 2 — 1.275. The solid line represents the fit 
VJI = -10, 256 cos(tf + 86°) km s" 1 . 



FIG. 4: Accumulation of the kick velocity extracted at r — 
75 M (top panel) and binary separation (bottom panel) as a 
function of time for the case L/M 2 — 1.375, 4>+ = 4>+ — 45° 
in Fig. El 




FIG. 3: Magnitude of the kick velocity (top panel) and % of 
energy radiated (bottom panel) as a function of initial orbital 
angular momentum for spin orientation <f>+ — cj>- — 45° . 



tion of $ = 9 — (f>, where d measures the angle between 
the initial spin and linear momentum vectors. The solid 
line represents the fit V\\ = -10, 239 cos(tf + 86°) km s" 1 . 
This is the same angular dependence found in quasi- 
circular orbits [161 ]. Furthermore, as with circular inspi- 
rals, the maximum kick is obtained for ■& rs 90° although 
in these hyperbolic encounters the kick is significantly 
larger, close to 10,000kms _1 . 

To investigate the dependence on the initial orbital an- 
gular momentum, we selected the case <f> + — 0_ = <fi = 
45° or -d = 108° in Fig. [2j which yields a kick magni- 
tude of 9, 589 km s~ and 15 % energy radiated, and car- 
ried out simulations varying the initial momentum of the 
BHs. Figure[3]shows the magnitude of the kick velocity V 
(top panel) and the energy radiated E^^ m % of the ini- 
tial energy (bottom panel) as a function of L/M 2 . Note 
that the maximum of V does not occur for maximum 
EVad- The case with the largest initial angular momen- 
tum, L/M 2 = 1.375, has an interesting feature as dis- 
played in the top panel of Fig. 5] There is a pronounced 
anti-kick before the recoil reaches its final value. The 
reason for this anti-kick could be due to the fact that the 
plunge is not as pronounced, appearing more circular- 
like. That is, there is a decrease in the rate at which 
the binary comes together, as one can see in the bottom 
panel of Fig. 2]during the time interval 20 — 40 M. Thus, 
the flux vectors responsible for the kick have more time 
to undergo the phase shift needed for the appearance of 
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FIG. 5: Magnitude of the kick velocity as a function of 
the extraction radius for different resolutions in the case 
(P/M = 0.2665, 4>+ =45°). From bottom to top are respec- 
tively resolutions of h = {M/0.4, M/0.5, M/0.63, M/0.8} at 
the extraction level. Dots are the data and lines correspond 
to the fit V(r) = Vi 2) + (Ki/r + K 2 /r 2 ). 



an anti-kick [22J. 

In addition to the errors from extracting the grav- 
itational recoil at a finite radius, the values of the 
kicks are affected by numerical finite differencing res- 
olution. To investigate these errors, we selected 
the case P/M = 0.2665,0+ = 45° from Table H 
and carried out simulations with resolutions h = 
{M/0.4, M/0.5, M/0.63, M/0.8} on the extraction level. 



Figure [5] shows the data (points) and the corresponding 

fitting (lines) to V{r) = v4? + K x /r + K 2 /r 2 . We find 
that at a given extraction radius, the error in the kick 
scales as ft 3 , and the kick itself grows with resolution. 
Based on this h 3 behavior, we estimate the error in the 
kick extrapolated as r — > oo computed using h = M/0.8 
to be of the order of a few hundred km s _ 1 , and we expect 
all the kicks presented to be of a similar accuracy. 

We have carried out a study of the gravitational recoil 
of the final BH in the merger of hyperbolic BBH encoun- 
ters. Wc have found that in general the kick velocities 
for in-plane initial BH spins arc significantly larger than 
those from the corresponding quasi-circular mergers. Our 
results suggest that kicks as large as 10, OOOkms^arc 
possible. We have also found that the dependence of the 
kick on the initial magnitude of the BH's spins is sim- 
ilar to the quasi-circular case. An analytic multipolar 
analysis of encounters for similar configurations can be 
found in Ref. [23||. A recent study by O'Leary, Kocsis 
and Loeb [28j has found that in dense population envi- 
ronments, there is a non-negligible probability for close 
fiybys or hyperbolic encounters [29||. Most of the cases 
they considered are those in which after the first pas- 
sage, the BHs release enough energy to become bound 
with large initial eccentricity. The hyperbolic encoun- 
ters we consider in the present work are the extreme case 
of immediate merger. Nonetheless, kicks with the magni- 
tudes found in the present study could lead to interesting 
astrophysical consequences. 
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